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We study the opening of the Haldane gap in two-leg and four-leg anisotropic spin ladders using 
bosonization and renormalization group methods, and we determine the phase diagram as a function 
of the interchain coupling and relative anisotropy, J z /J xy . It is found that the opening of the 
Haldane gap is qualitatively different for the two cases considered. For the two-leg ladder the Haldane 
gap opens for arbitrarily small interchain coupling, independent of J z /J xy , and the Haldane phase 
exists in a large region of parameter space. For the four-leg ladder the opening of the Haldane gap 
is strongly dependent on both the interchain coupling as well as J z /J xy , and the Haldane phase 
exists only in a narrow region about the isotropic antiferromagnet. 
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I. INTRODUCTION 

Recently the properties of systems which can be mod- 
eled by spin ladders have been extensively studied both 
experimentally and theoreticallyEl. It is well established 
that the magnetic properties of these materials depend 
strongly on the number of legs in the ladder. For even- 
leg ladders the susceptibility vanishes exponentially at 
low temperatures, while it shows power-law behavior for 
odd-leg ladders. This difference is due to the fact that 
the spectrum of magnetic excitations is gapless for odd- 
leg ladders and gapped for even-leg ladders. 

This resembles the alternance of a gapless and gapped 
spectrum for isotropic antiferromagnetkyspin-S' chains, 
where according to Haldane's conjecturecl the spectrum 
is gapless if S is a half-odd-integer while the spectrum 
is gappedj-for integer S. The two problems are not 
unrelatedBa, since a spin-S" chain can be described as 2S 
coupled spin-1/2 chains provided the interchain coupling 
is appropriately chosen. 

Although the appearance of the Haldane gap in an- 
isotrop>icr-|Spin-S' chains has been studied by several 
groupsa rW<pst earlier works on coupled chain or lad- 
der modelaljilil considered only the case where the cou- 
pling between the spins is isotropic. Some exceptions are 
the work by Watanabe et al&B using, bosonization tech- 
niques, and by Legeza and Sab/omll3 where the density 
matrix renormalization groupO was used to determine 
the phase diagram of anisotropic two-leg ladders. Since 
the results of the numerical calculations were not reli- 
able enough to get a definitive answer, in this work we 
use analytic methods to study the opening of the gap 
in a ladder model when two or four anisotropic spin-1/2 
chains are coupled with anisotropic interchain couplings. 
Our model is such that when the interchain coupling is 
equal to the coupling along the chains, the two-leg (four- 
leg) ladder behaves as a spin-1 (spin-2) chain. 

The rest of the paper is organized as follows. In Sec. II 
the Hamiltonian of the spin ladder model is presented, 
and its relationship to the composite spin representation 
of spin-S* chains is discussed. The model is bosonized in 
Sec. Ill, and the results known for spin-1/2 chains are re- 



capitulated. The two- and four-leg ladders are analyzed 
in Sees. IV and V, respectively, and the phase diagram 
as a function of anisotropy and interchain coupling is de- 
termined. Finally, in Sec. VI we summarize and present 
a discussion of our results. The paper also contains two 
Appendices. The first describes our bosonization conven- 
tions, while the second gives a derivation of the renormal- 
ization group equations. 



II. COUPLED SPIN CHAIN MODELS 



Consider a system of p (p = 2 or 4) coupled anisotropic 
spin-1/2 chains. Denote the spin operator on chain A 
(A = 1, ...p) at site j by Sj t \. The spins along the chains 
are coupled by a nearest neighbor Heisenberg exchange; 
the Hamiltonian for chain A is 



H x = j*yJ2[ s h s U,x 

+ jz ^2 s i\ s j+i,>. ■ 



(i) 



The coupling between chains can be taken in various 
forms. In spin ladder systems it is most natural to as- 
sume that the dominant interchain coupling acts between 
spins on the same rung and therefore has the form 
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8 h*h' 



(2) 



This kind of the interchain coupling is shown schemati- 
cally in Fig. [l| 

When the interchain coupling is strong enough and fer- 
romagnetic, the ground state and the low lying excited 
states of the ladder are those in which the p spins on rung 
j add up to form a larger spin of length p/2. Thus in 
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FIG. 1. Intra- and interchain couplings in a two- leg ladder 

this limit the two-leg ladder is equivalent to a spin-1 
chain in some of its physical properties, while the four-leg 
ladder gives an effective spin-2 chain. Strictly speaking 
this equivalence is valid only in the limit when the rung 
coupling becomes very strong. Nevertheless, it was es- 
tablished that in the isotropic two-leg ladder the Hal- 
dane gap-, is generated for arbitrarily small interchain 
couplingtij, i.e., its critical value is zero. 

Alternatively, one could n use the composite spin rep- 
resentation for spin chainsa. This representation is ob- 
tained by taking a single spin-S* chain with Hamiltonian 

H = ryY,(s*s* +1 + s?s? +1 ) + j*Y,s^Ui ( 3 ) 



and writing the spin operator as a sum of p = 2S spin-1/2 
operators 

S? = stl+sh + ... + s * p . ( 4 ) 

Inserting this into the equation above, we get 
p 



H — ^ H\ + ^ H\\' 

A = l A, A' 

A< A' 



(5) 



where H\ is the same as Eq. (g), so the first part is the 
Hamiltonian of p spin-1/2 Heisenberg chains, while the 
interchain coupling between them has the form 

H\\> = J± y ~^^ ( s j,A s j+l,A' + s |,A' s i+l,A 



S l\ S j+l,X + S j,A' S j + 1,A 



two spin-l/2's are added, in addition to the symmetric 
combinations that give the triplet 5 = 1 state, the anti- 
symmetric combination can also appear that corresponds 
to a singlet, S = state. It has been shownB, however, 
that the ground state and the low lying excited levels 
of this composite-spin model are identical to that of the 
true S — 1 model. The configurations where an S = ap- 
pears on at least one site have much higher energies than 
those we are interested in, since they effectively break the 
chains. Although it has not been checked numerically, it 
is reasonable that the same holds for the four-chain case, 
i.e., the lowest lying levels are those where the model 
is equivalent to the spin-2 chain; configurations where 
S = or S = 1 appear on at least one site effectively 
break the homogeneous chain, and therefore have higher 
energies. 

In actual ladder materials, in contrast to the cases dis- 
cussed above, the interchain coupling seems to be an- 
tiferromagnetic and act between nearest neighbor spins 
on the same rung. However, White showed that the two- 
leg ladder with (isotropic) antiferromagnetic coupling be- 
tween nearest neighbor pains is in the same universality 
class as the spin-1 chainlla. 

White's arguments were the following. Start with the 
composite spin representation for the spin-1 chain. Take 
"chain 2" and slide it by one lattice spacing to the left. 
Now, gradually turn off the coupling between sites y/E 
away; denote this coupling by J' ± . White was able to 
show that the system with J' ± = evolved continuously 
from the system with J' ± = J±. (i.e., There was no 
changein symmetry and no disappearance or appearance 
of gapsfea.) (See Fig. |.) 

White's arguments can be applied to the four-leg lad- 
der. For now, let us only consider isotropic couplings. 
(The case of anisotropic couplings will be discussed later.) 



(a) 



+ J 1 ^2 ( s i,A s j+i 



A' ' °j,X°j+l,X) 



(6) 




with 



j x y = j*y , ji = ,r- 



(7) 



One can view this, as shown schematically in Fig. as 
a ladder in which the interchain coupling is along diago- 
nals. 

jxy TZ 
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FIG. 2. Intra- and interchain couplings in a composite-spin 
model 

In this composite spin model the Hilbert space is larger 
than in the original spin-5 model. For example, when 



(b) 



(c) 



— s 

o e 


\— < 

5 — * c 


1 

> 


^ « 

i ^1 


» O 

) — ^-n 



T 6 

• e 


> 1 

> < 


► < 

► < 


> < 

> < 


) • 

3 O 



FIG. 3. White's argument showing that the (isotropic) an- 
tiferromagnetic two-leg ladder and the spin-1 chain are in the 
same universality class, (a) Start with the composite spin 
representation for a spin-1 chain, (b) Slide "chain-2" by one 
lattice spacing to the left, and gradually turn off the coupling 
between sites \fh away, (c) We are left with an antiferromag- 
netic two- leg ladder. 
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Start with the composite spin repersentation for the spin- 
2 chain. Shift "chain 2" and "chain 4" to the left by one 
lattice spacing. Then, gradually turn off the interchain 
couplings between sites which are not nearest neighbors. 
We are left with an antifcrromagnetic ladder with only 
nearest neighbor couplings. Therefore, as long as we have 
no change in symmetry and no disappearance or appear- 
ance of gaps as we turn off the couplings, the antiferro- 
magnetic four-leg ladder will be in the same universality 
class as the spin-2 chain. It seems reasonable that the 
two are related by continuity. Numerical studies have 
shown the antifcpmnagnetic four-leg ladder to have a 
gapped spectrumEBEy. It is hard to believe that, while 
performing the transformation above, we could go from 
one gapped phase to another independent gapped phase, 
and cross a critical point in the process. Therefore, it is 
reasonable that the four-leg (isotropic) antiferromagnetic 
ladder and the spin-2 chain are in the same universality 
class. 

One could ask how the Haldane phase of the ladder 
arises when the interchain coupling is switched on be- 
tween anisotropic spin- 1/2 chains which are gapless for 
— 1 < A < 1 and become antiferromagnetic with finite 
gap for A > 1. It is known that for anisotropic chains 
with integer spin, the Haldane phase shrinks to narrower 
and narrower ranges around the isotropic antiferromag- 
netic point as the spin length increases. Let A = J z / J xy . 
Denote the lower boundary of the Haldane phase for spin- 
S by A c i(S f ) and the upper boundary by A c2 {S). Then 



A cl (5i) < A cl (5 2 ) < 1 < A c2 (S 2 ) < A c2 (5i) 



(8) 



for Si < S 2 . Numerical calculations give A c i(5 = 1) = £ 
(e > 0, e < 1) and A cl (S = 2) sa 0.9. In the con- 
tinuum model of the spin-1 chain, e — is obtainedtiJ. 
The phase boundaries of the Haldane phase are shown 
schematically for the spin-1 and spin-2 chains in Fig. ^. 
Therefore, the interesting range of anisotropy for the lad- 
der is < A < 1, where the spin- 1/2 chain is gapless, 
while the spin-1 and spin-2 chains are in the Haldane 
phase for some range of A. 

Starting from the gapless situation and coupling the 
chains, taking into account that at the isotropic anti- 
ferromagnetic point the critical value of the interchain 
coupling is zerollHI, there are two natural ways 
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FIG. 4. Lower and upper critical anisotropies for the Hal- 
dane phase for S = 1 and S = 2. 
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FIG. 5. Possible scenarios for the opening of the Haldane 
gap. (a) & (b) Scenarios for the two-leg ladder, (c) & (d) 



Scenarios for the four-leg ladder. 
Haldane phase. 
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for the opening of the Haldane gap as the interchain cou- 
pling is switched on, as shown in Fig. ^|. Our aim is to 
decide which scenario occurs. For computational reasons, 
for the rest of the paper we will work with the composite 
spin model. In order to study how our system evolves 
from two (four) decoupled spin- 1/2 chains into a ladder 
that is equivalent to a spin-1 (spin-2) chain, as well the 
roles of the planar and Ising like parts of the interchain 
coupling individually, the Hamiltonian in Eq. (j^) will be 
generalized by allowing for ^ J xy and J]_ ^ J z in 
Eq. (|). 



III. BOSONIZING SPIN CHAINS 

While the spin-1/2 Hcisenberg chain can be solved ex- 
actly by the Bethe-Ansatz, models with higher spin val- 
ues can be treated, even in one-dimension, only approx- 
imately. One such approximation is to express the spin 
model in terms of fermions or bosons and develop a scal- 
ing theory for these degrees of freedom. In this section, 
we outline the key ingredients needed to fcrmionize and 
bosonize spin models. 



A. The Spin-1/2 Chain 

First consider a single anisotropic spin-1/2 chain de- 
scribed by the Hamiltonian (|l|), where the chain index A 
will be dropped. The Jordan- Wigner transformation 

s+ = sj + is) = c] cxp J2 c \ c i J ' 



s j - isV j = ex p (^ in J2 c l c i J s 



(9) 
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2 ' 



or because of the fermionic nature of the Cj , cj operators 
the equivalent form 



Sj = cjcos 




COS I 7T 2^ C]C; 
V 2=0 



(10) 

allows us to write this Hamiltonian in terms of spinlcss 
fcrmions 



+J' 



E 



(4 
(4 



C j+l C J 



C j~2 



(11) 



After Fourier transforming to momentum space 
H = J xy J2 cos ( ka ) c l c k 



J 



— J2 C0S (ia)p(q)p(-q) ■ 



(12) 



where a is the lattice constant, N the number of sites in 
the chain, and 



k 



/ j C k+q C k ■ 



(13) 



In the spin liquid phase where there is no net magnetic 
moment, the band is half filled. 

To study the low-energy properties of the model, we 
can linearize the spectrum around the Fermi points, 
fcj? = ±7r/2a. In this way a Luttinger model is obtained, 
where vf = J xy a is the unrenormalized velocity. In 
what follows we will work with a rescaled Hamiltonian, 
dividing H by J xy a, and use a dimcnsionless coupling 
A = J z /J xy (the relative anisotropy). 

Since the low-energy excitations of the Luttinger liq- 
uid are density fluctuations (which are bosons), it is con- 
venient to use, in the continuum limit, a further trans- 
formation from fermions to bosons. For clarity and com- 
pleteness, we have given a description of our bosonization 
conventions in Appendix A. 

For — 1 < A < 1 where the spin- 1/2 Heisenberg chain 
is known to be critical, the continuum bosonized form of 
the Hamiltonian is 



H = - I <Lr 



KW 



(14) 



For | A |<C 1 the parameters u and K can be determined 
perturbatively. However, the general form of u and K 
can be obtained by comparison with the Bethe-Ansatz 
solutior£3S. They are given by 



K 



2 (7r — arccos A) 



2 arccos A 



(15) 



An important property is that for A = (i.e., the free 
fermion point) K = 1, and for A = 1 (i.e., the isotroic 
point) K = 1/2. 

We will also need to know the bosonized form for the 
spin operators. Using the bosonization rules of Appendix 
A, in the continuum limit they can be shown to be 



s+(x) 

8>(Z) 



exp (— i^nQ) 



3 

a 



V7r it a 



e -i{-Kx/a) + C0S (V4^$) , 
(16) 



B. Coupled Spin-1/2 Chains 

For coupled spin-1/2 chains the same procedure can 
be used to express the spin operators in terms of boson 
fields, by simply attaching a A index to the field. The in- 
trachain part of the Hamiltonian will have the same form 
as Eq. (14). After rescaling by the bare velocity, J xy a, 
we write the Hamiltonian of the interchain coupling term 
in the form 



i v y 



S j,A' S j + 1,A 



/ + S 



j,X' a j 



+i,a) , (17) 



where we have let 

x <J xy and Jj - 



H 



A, A' 



jxy J- Jxy 

Using the bosonized form for the spin operators, we get 

dx r / , — , 
(2^i 5lCO H^^ + 

+ 5 2 COS (V4^($ A - $ A ') 

+.93 cos (\/^(e A -6v))] 
2JI 



+- 



-34 



+35 



dxd x <f>\d x $x> 
dx 



(18) 



(27TO) 2 

x cos (yiTr($\ + $ 



(v^(0a-6v)) 



A' 



dx 



cos (V^(e A - e v )) 



where 



93 



(2ira) 2 

X COS ^\/47T ($A - ^A 

3i = iJf , g 2 = -4J| , 
= -4tt J'^ , 34 = .95 = 2tt J x ^ . 



(19) 



Note that the coupling between the x and y compo- 
nents of the spins on different chains is non-local. 
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IV. THE TWO-LEG LADDER 

We first consider how the Haldane gap is generated in 
the anisotropic two- leg ladder when two gapless spin-1/2 
chains are coupled by an interchain coupling of the form 
given in Eq. (p) . The Hamiltonian of the system is 



H — Hi + H2 + Hi 2 



(20) 



As mentioned before for J^ y = 1 and J*[ = J z , our 
Hamiltonian is equivalent, as far as the low- lying states 
are concerned, to the true spin-1 Heisenberg chain. 
Therefore, we have a model which interpolates between 
two uncoupled spin-1/2 chains and a spin-1 chain. 
The bosonizcd form of the Hamiltonian is 



H= - I dx 



KUj+hd x ^i) 2 

A 



ATI:; 



+ - 



2J Z 



dxd x <$>id x <& 2 



dx 



gi cos 



(27ra) 2 
+ 92 cos (^/Att(^i - $2 
+ g 3 cos(V^(ei -e 2 ))] 



+ 



dx 



[34 cos (V7r(©i - ©2)) 



(27ra) 2 
x cos ^V47r(<i>i + $2) 
+ g 5 cos (V^(Bi - 2 )) 
x cos ^\/47r(<i>i - $2) 

It is useful to define the fields 



(21) 



1 



($1 + $2) , $0 



1 



($1 - $ 2 



In terms of these fields our Hamiltonian has the form 



(22) 



H = y I </-<■ 



A" S H 2 + -kdA) 2 

As 



-01 



dx 



(277a) 2 



dx 



K a U 2 a + —(dA) 2 



52 cos v 87r$ a + .g 3 cos 



(27ra) 2 

~ ff4 / (2^ C ° S ( V ^ 6a ) C ° S (^ s ) 
-.95 



where 



A" s = A- 1 + 



2KJ Z 



K„ = K 1 



2a"J* 



-1/2 



-1/2 



« s = It 1 + 



2a" J? 



1/2 



U a = U 1 



2a" Jf 



(24) 
1/2 



For J x ^ ,J Z ^ Clwe have 



A" s « K \1 - 
K n *k(i 



KJ Z 



UTT 

KJ Z ± 

UTT 



, u s pa u ^1 + 

, «„ W U f 1 - 



A"J1 
U7T 

UTT 



(25) 



y (2™) 2 ^° <v^v vu "~v 

/(2^ C0S ( V ^ ea ) C0S ( V ^ <i>a ) ' (23) ^ Z 



We are interested in whether or not the interchain cou- 
pling causes a gap in the excitation spectrum. Therefore, 
we would like to identify the relevant operators; these 
operators will "pin" their arguments, thus causing gaps 
to appear. To do this we consider the scaling dimen- 
sions of the operators in the interchain coupling (see Ap- 
pendix A). The scaling dimension of cos (\/8tt<& s ) is 2a" s ; 
the scaling dimension of cos (\/87r$a) is 2A" a ; the scal- 
ing dimension of cos (-\Z27r0a) is l/(2A" a ); the scaling di- 
mension of cos (V27r9 Q ) cos (\/4ir§ s ) is 2a" s + l/(2K a ); 
the scaling dimension of cos (y27rQ a ) cos (\/47r$ a ) is 
2A" a + l/(2A" a ); Therefore, gi will grow at large distances 
for K s < 1; g 2 will grow for K a < 1; g 3 will grow for 
K a > 1/4. The f/5 operator is always irrelevant, and 
therefore will be ignored. However, we will keep the 34 
term since, though seemingly irrelevant, it is the most 
relevant operator generated by the interchain coupling 
which couples the symmetric and antisymmetric modes. 
We will see that the 54 term plays a subtle role. 

It is interesting to study the behavior of the xy and z 
components of the interchain coupling individually. First 
we consider the case J x ^ — and JJ ^ 0. Therefore, 
only gi and gi are nonzero. The gi term is relevant for 
< A < 1. The <7 2 term is relevant for K < 1 — — 
(for weak interchain coupling.) Therefore, the symmetric 
mode is gapped in the entire region < A < 1, while the 
antisymmetric mode remains gapless in a region about 
A = and is gapped outside of this region. In the re- 
gion where $ Q is-pinned, the antisymmetric mode is in 
an ordered phaseH. 

Next we consider the case J 3 ^ ^ and J Z L =Q. There- 
fore, only g 3 and 54 are nonzero. For < A < 1, the 53 
term is relevant while the 174 term is irrelevant. There- 
fore, 53 will grow while 54 will initially decrease under 
the RG. However, once gs grows to 0(1), we can safely 
say that Q a is pinned. At this point, assuming that g^ is 
renormalized to 54, the term that couples the symmetric 
and antisymmetric modes can be written as 



dx 



(27ra) 2 
Jr { cos ( \ 2^(-)., ) ^ 



(V^Ga) ^ COS (V87T$ a ) 

(V8^$ s ) . (26) 



dx 
{2na) 2 
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TABLE I. Results for the two-leg ladder in the region < A < 1 (i.e. 1 > K > 1/2). 



J x v = 0, Jl £ 


Jl y ^ 0, J| = 


JI/Jf = A 


<£> s pinned 


pinned 


pinned 


$a, ©a $ a pinned for K a < 1, 


8 a pinned, 


a pinned, 


ordered 


disordered 


disordered 


Gapless for X a < 1, 






Phase Haldane for X a < 1 


Haldane phase 


Haldane phase 



The expectation value is taken with respect to 

1 



Ua 

2 

IB 



dx 
dx 



{2-Kaf 



cos ( v2vr6a 



(27) 



where K a and 173 are the renormalized values of K a and 
173. Then we can define a new coupling constant 



<?4 



94 { cos 



(28) 



Therefore, our effective Hamiltonian for the symmetric 
mode is 



H cS = — I dx 



—2 1 



K s ui + ={d x <s> s y 



94 



dx 



(27ra) 2 



8tt$. 



(29) 



where K s is the renormalized value of K s . Now 
cos (v / 87r < i ) s) is relevant for K s < 1. Therefore, the term 
is relevant. We see that J xy alone generates a gap in both 
the symmetric and antisymmetric modes. However, since 
Q was. pinned, the antisymmetric mode is in a disordered 
phase B 



Finally, we consider the case ^ and J]_ 7^ 
with J^/J^ — A. Therefore, for J x ^ — 1, we recover 
the results for the anisotropic spin-1 chain. The g\ and 
g 3 terms are relevant in the entire range < A < 1, 
and 53 is always more relevant than gi- We see that 
both the symmetric and antisymmetric modes will get 
gapped, with the antisymmetric mode being in a disor- 
dered phase; the Haldane gap is generated for arbitrarily 
small interchain coupling. 

The results for the two-leg ladder are summarized in 
Table |. 



V. THE FOUR-LEG LADDER 

Next we study the opening of the Haldane gap in the 
four-leg ladder. The Hamiltonian is 

H = H 1 +H 2 + H 3 + H 4 

+Hl,2 + #1,3 + fll,4 + ^2,3 + H2.4 + H3A ■ (30) 

Let us introduce 3> T = ($1, $ 2 , $3, $4) and n T = 
(Hi, n2, n3, n4). Then in its bosonized form our Hamil- 
tonian is 



H = j dx 
+.9i , 



^n T n + (d x *) T M(d x * 

r (2^H^ ($2+<f>i 



+ .92 



93 



+ 94 



cos (yjir(^ 3 + $1)^) + cos ^V47t($ 4 + $1) ] 
h cos (yin(^ 3 + $ 2 )) + cos ^V47r($4 + $ 2 )) + cos (\/4tt($ 4 + $ 3 )) 
7^x2 cos (^4tt( $ 2 - $1)) + cos (%/4tt(<I>3 - $1)) + cos (%/47t($ 4 - *i) 
f- cos ^V47t(<I>3 - <5> 2 )) + cos ^\/47r(<I>4 - $ 2 )) + cos ^V47r($ 4 - $3) 

' / (2to) 2 ^ C ° S ^^ 02 ~ + C ° S (^( 9 3 - 9l )) + C0S (^( 9 4 - 

+ cos (V^(e 3 - e 2 )) + cos (V^(e 4 - e 2 )) + cos (V4^(e 4 - e 3 )) 

t J , 2 ^ a y C0S (V^(02 - ©l)) COS (%/4^($ 2 + + COS (y/n(9 3 - ©l)) COS (V4tt($3 + 

cos (Vt?(©4 - ©l)) cos (v / 4tt($4 + $1)) + cos (Vt?(©3 - ©2)) cos (V4^($ 3 + $ 2 )) 



©l)) 



G 



+ 95 



+ COS 

f dx 



(V¥(e 4 - e 2 )) cos (V44$4 + $2)) + cos (VtF(B4 - e 3 )) cos (y4*(® 4 + $3)) 

(V^(©2 - 0l)) cos (V4^($ 2 - + cos (V^(9 3 - 9l)) cos (V4^($ 3 - $x)) 

-I- ens 



(2nay L " V V 

+ COS (V7r(6 4 - 6l)) COS ^V47r($4 - $1)) + COS (V7r(9 3 - 6 2 )) COS ^\/47r($3 - $ 

- (^(©4 - 62)) cos (V4^($ 4 - $2)) + cos (V^(0 4 - 63)) cos (\74^($ 4 - $3)) 



)) 



+ COS 



(31) 



where M is given by 



M = 



/ a b b b\ 

b a b b 

b b a b 

\b b b a J 



(32) 



with a = u/(2K) and b = J]_/tt. M is diagonalized by 
the orthogonal matrix S 



(33) 



It should be noted that the choice of S is not unique. 
Define A by 





fj 


1 

1 


1 
1 


*\ 


s = 


I 


I 


I 


I 




u 


1 

2 


1 
2 


-I J 



A = S J MS. 



A is given by 



/ a + 'ib 
a- b 
a-b 



\ 





V 0-6/ 
Define the fields $ a i, $ a2 , and <J> a3 by 

/I I I I \ 
/ ? 2 2 2 \ 



(34) 



(35) 



1 1 

? ! _! _? 

I _! _f I 

I _! ? J 
2222 




(36) 



In terms of the old fields, the new fields are given by 



<S> S = I ($1 + $2 + *3 + $4) 

= I ($1 + $2 - $3 - $4) 

$a2 = I ($1 - $2 - $3 + $4) 

*a3 = I ($1 - $2 + $3 - $4) • 



(37) 



We see that $ s is a symmetric mode while $ a i, $ a2 , and 
$ a 3 are antisymmetric modes. 

Let us define * T = ($ s , $ ai , $ a2 , <& a 3) and fl T = 
(n s ,n ol ,n o2 ,n a3 ). Then $ = S* and n = sri. In 
terms of the new fields, the quadratic part of the Hamil- 
tonian is given by 



dx 



(38) 



We can write this as 



H a = — I dx 



K s nl + — (d x <S> s f 



3 . 



KaUl l + — {d x $ ai f 



(39) 



with 



U s = U 1 + 



UTT 



U a = U 1 - 



2ifJf 



1/2 



1/2 



if s = if 1 + 



if a = if U 



6if 


r. 








(40) 


2if J| N 


r- 


U7T , 





For J^ y , J| « 1 we have 



u. w u H 



««*» 1 



UTT 

KJ 



UTT 



X V if„«if (l + 



3if Jj 
if Jj; 

UTT 



(41) 



In terms of these new fields, the interchain coupling is 
given by 



'1 J ^a) 2 C ° S (^" $al ) cos (^4tt $ . 

(\/47r$ a2 ) cos (v^^al 



2ff 

+ 2 52 / 



dx 
a) 
dx 

dx 



cos 
cos 



^V47T$ a2 ^ cos (VJn^^ + cos ^V47r$ a3 ^ cos (\/47r$ 



(V47r$ a3 ) cos (V4K<i>ai) + cos (V47r$ a3 ) cos (\/47r$ a2 ) 

+2.f/3 / ^n)2 [ cos (^©02) cos (VnO a i) + cos (V7re a3 ) cos (Vn@ al ) + cos (VwQca) cos (Vtt©. 
dx r / .— ... ^ x\ / at— I ,\ / r- 



J o2 )] 

: 34 y 7^)2 C0S (V^(0a3 + ©a 2 )) COS (\/47($ s + $ al )) + COS (V^(6 a3 - 6 a2 )) COS (V^($ s - $ Q i)) 
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+ cos (V¥(e a2 + o i)) cos (VI^(<s> s + $ o3 )) + cos (v^(e a2 - e i)) cos {y^{<s> s - $ a3 )) 

+ COS (V^(6a3 + G Q i)) COS (VI^(<S> S + $ o2 )) + COS (y/w(Q a3 - o l)) COS (VA^($> S - *a2)) j 

cos (V^( «3 + ©02)) cos (yi^(<P a3 + $ o2 )) + cos (V^(G a3 - G a2 )) cos (\/47r($ o3 - $ o2 )) 

1 COS (V¥(0 a2 + G Q i)) COS (V44$ a2 + $ a l)) + COS (V^(G a2 - G a l)) COS (\/I^($a2 - $al)) 

cos (v^r(Q a i + G a3 )) cos ^V4^($ a i + $ a3 )) + cos (V¥(G a3 - G a i)) cos (\/4^($ a3 - $ ai ) X 



+ COS 

+ 



(42) 



We see that gi and #4 couples the symmetric mode to the 
antisymmetic modes, while g 2 , g 3 , and 175 only couple the 
antisymmetric modes among themselves. 

Just like for the two-leg ladder, we would like to know 
the relevance or irrelevance of the operators in the in- 
terchain coupling. Therefore, we consider the scaling di- 
mension of the operators in the interchain coupling. The 
gi terms have scaling dimension K s + K a ; the gi terms 
have scaling dimension 2K a ; the 33 terms have scaling 
dimension l/(2K a ); the g± terms have scaling dimension 
K s + K a + l/(2K a ); the #5 terms have scaling dimension 
2K a + l/(2K a ). Therefore, the g\ terms are relevant for 
K s + K a < 2; the gi terms are relevant for K a < 1; the 
53 terms are relevant for K a > 1/4. The g§ terms are 
always irrelevant and will be dropped. However, we will 
keep the #4 terms. Though they are irrelevant, they are 
the most relevant terms (arising from the xy part of the 
interchain coupling) which couple the symmetric and an- 
tisymmetric modes. We will see that the 54 terms play a 
subtle role. 

Since the g\, gi, and 53 terms are all relevant and all 
of the modes are coupled together, the physics is deter- 
mined by the operators which reach strong coupling first 
under the RG. The RG equations for the parameters are 

^f- = [2 -(K, +*„)] Si 



dg3 
dl 

dgi 
dl 



1 
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2K a 
2-(K s +K a + 



2K n 



dK s 
dl 

dKa 

dl 



= -2.9? 



4ttu s 

( Kg 

\AnUa 



%2 2 



K s 

AlTUs 

Kg 

A-KUg 



+ 53 



Airug 



-A 1 K ° 



\Airu a 



2 ! / 1 



(43) 



.94 



See Appendix B for a derivation of the RG equations. 

Just like for the two-leg ladder, it is interesting to study 
the behavior of the xy and z components of the interchain 
coupling individually. First we consider the case J xy = 
and J\_ 7^ 0. For this case the initial values of the cou- 
pling constants are gi(l = 0) = 4Jf_, g 2 (l = 0) = — 4Jj% 
gs{l = 0) = 0, and g±(l = 0) = 0. The g\ terms are 
relevant for < A < 1. The g 2 terms are relevant for 
K < 1 — — (for weak interchain coupling.) However, 
since all of the modes are coupled (and hence locked to- 
gether), both the symmetric and antisymmetric modes 
are gapped in the entire region < A < 1. Notice that 
since the $ a i were pinned, the antisymmetric modes are 
in an ordered phase. 

Next we consider the case J x ^ 7^ and J 2 ^ = 0. For 
this case the initial values of the coupling constants are 
9l (l = 0) = 0, g 2 (l = 0) = 0, g 3 (l = 0) = -AirJ 3 ?, 
g A {l = 0) = 2nJl y , K s (l = 0) = K, K a (l = 0) = K, 
u s = u, and u a — u. The g 3 terms are relevant and 
the gi terms are irrelevant. Therefore, g 3 will grow and 
94 will initially decrease under the RG. However, once 
g 3 = O(l) we can safely say that the Q a i are pinned. 
Then we can write the g^ term as 



.94 



/ 



{2na) 



f (cos (V¥(G a2 + o3 ))> cos (y^(<P s + *oi)) + (cos (v^r(9 i + a2 ))) cos (V5r(* a + $ a3 )) 
+ (cos (Vtt (G a i + 6 a3 ))> cos (V4T($ S + $ a2 )) + (cos (V^F (G a3 - Gai))) cos (V4^($ s - $ a2 )) 
+ (cos (Vtt (G a2 - Gai))) cos (V4^(* s - $ a3 )) + (cos (Vtt (G a2 - G a3 ))) cos (\/4^($ s - $ Q i)) 



gl(cos(V^(G a i +e o2 ))> J 



COS ^V4^($ s + $al)) + COS (Vl^s + $ a3 )) 
+ COS (V44$ s + $ o2 )) + COS (yi^(<S> s ~ $02)) COS (V4n(® a ~ $a 3 )) + COS (V4n(® s - $ a l)) , (44) 

where ~gi is the renormalized coupling. The expectation value is taken with respect to 
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—„2 1 



^ 2?ro j 2 [cos (^©02) cos (VttQoi) + cos (V7re a3 ) cos (V7r@ai) + cos (VirO a3 ) cos (v/tt^)] , (45) 

where K a and 53 are the renormalized values of K a and 53. Notice that H a is invariant under Q a i — > —O a i] therefore, 
(cos (^/n(Q ai — @ a j))) = (cos (y/n(<d ai + Q a j))) ■ Now we can define a new effective coupling 



and write the 34 term as 



.9 4 = .94 (cos (Vn(Q al + ©02))) 



(46) 



(47) 



Recall that [Q(x), 0(y)] = 0(1). Since the <d a i are pinned, the 4> a i will fluctuate wildly and e 1 ' 3 *" will have 
exponentially decaying correlations. Looking at the interchain coupling, one might think that the g' 4 term is now 
irrelevant due to the presence of the cos(v / 4~7r < £'ai) terms. However, this is not entirely correct. To see what can 
happen we evaluate the (normalized) partition function perturbatively in g\. We write 



Z 



V^e~ So e~ s 



l-(S 1 ) + -(S 2 1 )a + 



where 



So = -=^= I d 2 x 
2K, 1 



— — I d x 



\{d T Q ai ) 2 + {d x Q ai f 



(27ra) 2 

with K s the renormalized value of K Sl ( )o denotes averaging with respect to Sq, and 

3 * rf. 



/d 
-fiAp COS (V4^$ s ) COS (y^t^ai) 
I 



The first non- vanishing correction is 



(s() - nmy 



2 f d 2 xi d 2 X2 



(27ra) 2 (2na) 



cos ^v / 47r$ s (a:i)^ cos ^V47r$ s (iz; 2 )^ cos ^V47r$ a i (xifj cos (^V4^$ a i(x2 



(48) 



f d 2 x 

2ffi J ^ 27ra )2 [ C ° S (^© a2 ) C0S (V^Qal) + C0S (Vn@a3) COS (y/nQal) + COS (^©03) COS (v / 7T©a2)] (49) 



(50) 



(51) 



(We have ignored correlations in fluctuations between 
<& a i and & a j for i ^ j, since they are suppressed.) 
Since e^* Qi have exponentially decaying correlations, we 
approximates 



^COS (yi77^ a i(xi)j COS {y^TT<$>al(x 2 )^ / 

Ri D (5(a;i — ir 2 ) + higher order corrections , 
where D is an 0(1) constant. Therefore, we have 

d 2 Xi 



(52) 



4 = i + §<V.>'d/ 



(27ra) 2 



cos' 



= 1 + 3(g' 4 f D | ■ 



d?Xi 

Re-exponentiating, we have 



cos 2v47r$ 



Z 



- < e - S; )o 



where 



d 2 x 

(27TO)' 



cos(Vl67r$ s ) . 



(53) 



(54) 



(55) 
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TABLE II. Results for the four-leg 


ladder in the region < A < 1 (i.e. 1 > if > 1/2). 




Jj_ — U, J j_ f= u 


T xy n J z — o 






pinned for < A < 1 


pinned for if s < 1/2 


pinned for if s < 1/2 




$ a pinned, 


© a pinned, 


8 a pinned, 




ordered 


disordered 


disordered 






Gapless for TT S > 1/2, 


Gapless for ~K~ S > 1/2, 


Phase 


Haldanc for < A < 1 


Haldane for K s < 1/2 


Haldane for ~K~ S < 1/2 



Therefore, our effective Hamiltonian for <J> S is 



H e ff = — I dx 



H e ff = -zr I dx 



K S U 2 S + ={d x § s y 

ft* 




(56) 



where we have introduced a new effective coupling con- 
stant, g. 

Now we have an effective Hamiltonian for <i> s which is a 
standard sine-Gordon Hamiltonian. The scaling dimen- 
sion of cos (v / 16~7r<I> s ) is 4K S ; it is relevant for K s < 1/2. 
Therefore, g grows strong for K s < 1/2; $ s gets pinned 
(and hence becomes massive) for K s < 1/2. However, the 
Q a i are pinned for < A < 1. Therefore, the antisym- 
metric modes are in a disordered phase for < A < 1. 

Finally, we consider the case J x ^ ^ and JJ ^ 
with J^/J^ y = A. Therefore, for J^ y = 1, we re- 
cover the results for the anisotropic spin-2 chain. For 
this case the initial values of the coupling constants are 
g x {l = 0) = 4JI, g 2 (I = 0) = -4J1, g 3 (l = 0) = -47rJf , 
and g±(l = 0) = 2kJ x ^ ' . Looking at the RG equations, 
we see that for 1 > K > 1/2 (i.e. for < A < 1), the 
gs term will dominate; it is the most relevant operator. 
Therefore, the Q a i fields will get pinnned. Once the Q a i 
fields are pinned, just like the J x ^ ^ & J\ = case, 
the (74 term becomes 



'Mi 



3 
i=l 



(27ra) 2 



. (V^$ 8 ) cos (%/4^r$ ai ) , (57) 



where g\ = 54 (cos (i/tt(@i + 62))), and ~gi is the renor- 
malized value of 34. Define g[ — 2(g~ 1 + g' 4 ). ('g 1 is the 
renormalized value of g\.) Then our interchain coupling 
is 

^interchain =9^ J ^a) 2 ^ (^^ $s ) C0S (v^^ai) 

(58) 

From this point, the analysis proceeds identically to 
the case ^ & J\_ — — derive an effective sine- 
Gordon Hamiltonian for the symmetric mode and deter- 
mine the phase boundary from this effective Hamiltonian. 
Similar to the case J xy ^ h, J]_ = 0, the effective sine- 
Gordon Hamiltonian is 



K S U 2 S + ={d x <5> s ) 2 



dx 



(27ra) 2 



(59) 



where K s is the renormalized value of K s , and g is an 
effective coupling constant. Again, $ s will become mas- 
sive for K s < 1/2, and the Q a i are pinned for < A < 1. 
Therefore, a Haldane gap appears for K s < 1/2 and 
the antisymmetric modes are in a disordered phase for 
< A < 1. 

The above results for the four-leg ladder are summa- 
rized in Table [n]. 

To determine the phase boundary between the gapless 
and Haldane phases for the case J x ^ ^ 0, J]_ ^ 0, with 
J^/Jj^ = A, we numerically integrated up the RG equa- 
tions. It turned out that g 3 grew so much faster than g\ 
or <72- Our numerical procedure was the following. We 
integrated up the RG equations to a scale, £. £ was de- 
fined as the scale where g 3 = 0(1) (g 3 ^ gi, g 2 , 34-) At 
that point, we could safely say that the Q a i fields are 
pinned (up to gapped fluctuations)E3. Then we defined 



(60) 



From the above arguments, K s — 1/2 defines the phase 
boundary. 

Unfortunately, we were only able to determine the 
phase boundary qualitatively. This is because K s is de- 
fined when g 3 — 0(1). (We need 33 = 0(1) so that we 
can safely say that the 6 Q i fields are pinned.) However, 
our RG equations are no longer valid when any of the 
couplings are of 0(1), since they were calculated to low- 
est order for weak coupling. Therefore, our procedure 
becomes uncontrolled. However, as long as there is no 
intermediate fixed point, our results will be qualitatively 
correct. Although we cannot think of what such a fixed 
point would physically correspond to, it cannot be ruled 
out. The phase diagram is shown in Fig. ^| 
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HALDANE PHASE 



GAPLESS PHASE 



FIG. 6. Phase diagram of the four-leg ladder. The solid line 
shows the phase boundary between the gapless and Haldane 
phases calculated from the RG equations. The dashed line 
indicates a speculated phase boundary in the region where 
the perturbative treatment is not valid. 



in Fig. g. Therefore, for the two- leg ladder it is scenario 
(a) of Fig. |^ which occurs, while for the four-leg ladder 
it is scenario (d) of Fig. ||| which occurs; the opening of 
the Haldane gap is qualitatively different for the two-leg 
and four- leg ladders. 

How do our results apply to spin ladders where the cou- 
pling is antiferromagnetic and between nearest neighbor 
spins on the same rung (as in actual ladder materials)? 
In Sec. II. we argued that the isotropic composite spin 
model has the same phase diagram as the isotropic anti- 
ferromagnetic ladder. Now what if we include anisotropy 
and allow for the interchain coupling to be different from 
the coupling along the chains? It is likely that the phase 
boundaries (between the massless and Haldane phases) 
in the A — Jj_ plane, which we found using the composite 
spin model, will shift. For the two- leg ladder it appears 
that the Haldane phase still opens from A = 0, J± = 0c3. 
However, for the four-leg ladder, it is possible that the 
Haldane phase only exists at the isotropic point or in an 
extremely narrow sliver about the isotropic point. We 
leave this (and other possibilities) for future work. 



VI. CONCLUDING REMARKS 

In this paper we studied the generation of the Haldane 
gap in two-leg and four-leg anisotropic spin ladders us- 
ing a particular form for the interchain coupling, namely 
the case where spins situated diagonally interact. An in- 
terchain coupling of this form was chosen so that when 
Jj_ = 1, the two- leg (four- leg) ladder was equivalent to a 
spin-1 (spin-2) chain (as far as their low lying excitations 
are concerned.) 

After bosonizing our model, we made a (linear) trans- 
formation to a new set of fields. One field contained a 
symmetric combination of the operators on the legs; the 
others were, in a general sense, antisymmetric combina- 
tions. For the antisymmetric modes, whenever a or G a j 
were present, the operators containing <d a or Q a i were 
the most relevant operators; the a or a j were pinned 
in the entire range < A < 1. However, when only $ a or 
Q a i were present, the behavior depended on the number 
of legs of the ladder. For the two-leg ladder there was a 
region about A = where <E> a remained gapless; for the 
four-leg ladder, the $ a j were pinned in the entire range 
< A < 1. 

For the symmetric mode, it was always $ s (rather 
than S ) which came into play. However, its behavior 
was strongly dependent on the number of legs of the lad- 
der. For the two-leg ladder the symmetric field, was 
pinned in the entire range < A < 1, and the Haldane 
gap was generated for arbitrarily small interchain cou- 
pling. However, for the four-leg ladder, the Haldane gap 
was strongly dependent on both the interchain coupling 
and A. In order to determine the phase diagram, the 
RG equations for the parameters were derived and inte- 
grated. It was found that the Haldane phase appeared 
in a narrow range around the isotropic point, as shown 
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APPENDIX A: BOSONIZATION DICTIONARY 

The bosonization procedure is well documentedElEi 
However, since various conventions exist in the literature, 
we present our conventions here for clarity and complete- 
ness. 

To study the low energy properties of our fermion 
model, we linerarize the dispersion about the Fermi 
points and write 

r(>(x) = c j /y/E = e-* k r^ R (x) + e ikFX i, L {x), (Al) 

where x = ja. .— . 
Then we have the following bosonization rulesa 




exp(iv47r</>_R) , 
cxp(-iv / 47T</>L) , 



ipli(x)tpR(x) + ^l(x)^ L {x) = po + -=d x $(x) 



(A2) 
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where the chiral boson fields are related to the standard 
Bose-field, and its dual field, 0, via 

$ = fa + (f>L , <3 = <j>R~(f)L- (A3) 

The fields $ and satisfy the commutation relations 

[$(x),e(y)]=i6(y-x), (A4) 

where 9 is the step-function. $ and can be thought 
of as order and disorder fielda n respectively, in the usual 
statistical mechanical sense.clUa 

Our fields are normalized such that the free boson ac- 
tion is given by 



So -2 



(d T $y + (^$) 



The corresponding Green's function is 



1 / R 2 



(A5) 



(A6) 



where R is an infrared cutoff (i? — > oo). It is the solution 
to 

- (d 2 T + d 2 x ) G(x, t) = 6(x - x )S(t - r ) . (A7) 

Note that in our calculations we use the lattice spacing, 
a, as an ultraviolet regulator. 
For a Hamiltonian of the form 



H 



dx 



KIT 



1 

K 



we make a canonical transformation 
1 



n = 



K 



n, $ 



(A8) 



(A9) 



In terms of these new variables, our Hamiltonian has the 
form 



H 



dx 



n 2 



(A10) 



Using this form, we can pass to a Lagrangian description 
and obtain an action of the form in Eq. ( |Aq ) . 

In the critical regime of the spin-1/2 Heisenberg chain, 
the spin-spin correlations functions have power law be- 
havior which can be calculated in the boson representa- 
tion using the relationsEj 



e ia$(x) e -ia<i>(y) \ _ 



\x - y\" 



where we have fi rst used the transformation in Eq. ( |AS| ). 
From Eq. (All) we can obtain the scaling dimensions of 
the operators e ia *W and e l/3e ^; it follows that e ia ^ 
has scaling dimension and e z ^ e ( x ) has scaling dimen- 

sion . An operator is relevant if its scaling dimension 
is less than 2; it is irrelevant if its scaling dimension is 
greater than 2. If an operator is relevant, its coefficient 
grows strong at large distances; if it is irrelevant, its co- 
efficient becomes smaller at large distances. Therefore, 
e ia<s>(x) - 1S re ;[ evan t f or °lJL < 2; e l 0®( x ) is relevant for 

< 2 



APPENDIX B: RENORMALIZATION GROUP 
ANALYSIS 



For weak interchain coupling we can analyze the 
physics of our model with a perturbative RG analysis. In 
an RG treatment we coarse grain our system, integrating 
out short wavelength (high energy) degrees of freedom, 
and obtain an effective theory for the long wavelength 
(low energy) degrees of freedom. 

We will derive the RG equa tions for <?i, g± and K s in 
detail to illustrate the mcthodEZla . The RG equations 
for the other parameters can be derived in a similar way. 

We consider the correlator (^v^M^e-iv^M^)) 
and evaluate it perturbatively in the interaction. There- 
fore, we write 



giV47r<I> s (xi)g— iV4'7r<I> a (x2) 



^iV^TT® s (X l) £ — iV4TT$s(x2) (J \ 

/ ,,t\/35r*,(xi) -iV37r*,(x 2 ) o2\ , 

~2Y 6 bl /o + "' 



(Bl) 



where ( )o denotes averaging with respect to the free bo- 
son action. The first nonvanishing correction comes from 
the Sf term. 



o i0Q(x) -i/30(y)\ = 



/ \\x-y\' 



(AH) 



/ „i\ZiTT^ s (xi) ~i\fiTT^ s (x2) q2\ 

V 6 bl /o 



d 2 X3 d 2 Xi 



i\/47r$ s (ii) — iv47r<I> s (x 2 ) ivi^^ s{x 3 ) — i\/47r<I> s (2:4) \ / iv47r$ a i (X3) — iV47r5> a i (x 4 ) 
o c o o / \ u o 



+ ~ 2 9t 



{2ira) 2 (2na) 2 

d X3 d X^ / p i v ^$ 3 (^ 1 ) p _i v ^$ s ( a ; 2 ) (= i V q7$ s ( 2 ;3) (= _i v ^7$ s ( 2 ; 4 )\ 1 ^ ^ j V^TT* a 1 ( X 3 ) p - 1 \TAtT $ a 1 ( X4. ) 

(27ra) 2 (27ra) 2 
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Evaluate the correlators to get 



(2ffi) 2 



d 2 X3 d 2 xn 
(2ira) 2 (27ra) 2 \J *2l I 



d 2 X3 d 2 X4 

(27ra) 2 (27ra) 2 \ \ z 21 \ 



" ^31 | 2 | Z 4 2 r K 
Z 3 2 | 2 | *41 I 2 



231 I 2 1 Z42 12 ' '"" 



a 



Z32 | 2 | ^41 | 2 



243 r 

2 x K 



K, 



I I 2 



K„ 



Z43 



where Zij — u s (ri — r,-) + — Xj) and z°- = u a (ri — tj) + i(xi — Xj). Define 



r, = u s Ti x + Xi y 



and make the change of variables 



I I 2 



(2) 



(3) 



(4) 



r = r 4 - r 3 



R = 2 ( r 4 + r -3) 



Then we have 



Z21 



A'., 



d 2 i? d 2 r 
(2^)2 (2W) 2 

K B +K a ( 



3 2 1 



^u 2 VI 221 



A'» 



d 2 i? d 2 r 
(27ra) 2 (27ra) 2 




-, A s 



R-±r-r 2 | 2 |R+ir- ri 



2X A s +A a+5 i- 



A„ 



sin 2 9 



cos 2 



where u s (t4 — T3) = r cos# and £4 — 2:3 = r sin 9. Write 



1 



V 



sin 2 0+ M cos 









V s 






= l + /3 


1 - 











cos 2 H . 



(5) 



(6) 



(7) 



Since 1 — (^^j — ®(J±)i the anisotropy between u a and u s gives us corrections which are of higher order. Therefore, 

to this order we ignore the anisotropy. 

Using that the most singular part of the integrand is for r near 0, we expand the integrand about r = 0. We get 



R-^r-n | 2 |R+±r-r 2 | 2 



R-±r-r 2 | 2 |R+i] 



K, 



= -4K 2 r 2 cos 2 * ( R - r i)-( R -^) 



where we have kept only connected piecesBS 
Therefore we have 



R-n | 2 | R-r 2 



disconnected pieces , 



(8) 



i\/47r& s (x) — i\/47r$ s (y) q2 



221 



A, 



ri 2 i? 

(27ra) J 



V R ln 



3 2 
2 54 



R 



•VrIii 



R r 2 



2iff 



(27ra) : 



-r 2 cos 2 ( 



-3ffi 



A s +A a + T 



2^; 



d 2 r 
(27ra) ; 



■r 2 cos 2 



A s +A a 



(9) 
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where we have used that 



< R -> =v R ,„f |R - r 



R-r 



Integrate over 9, integrate by parts over R and ignore the surface term, and use that 

R r 



V^ln 



= 27r<P(R-r) 



to get 



^21 



K„ 



In 



r 2 - >"i 



+ 3<?I 



65? 



dr /r \3-2(K s +K a ) 



dr /r\ 
a \a/ 



dr /rx 3-2(/f 3 +K a+T i-) 



dr /r\ 
a \aJ 



(10) 



Finally, we have 



^21 



A', 



1 + In 



r 2 ri 



3.9, 



j 

/dr /r\ 
a \a/ 



3-2(K 3 +K a ) 



Attu. 

dr /r\3-2(K s +K a +^-) 



/dr /ry 
a \a/ 



(11) 



Treat the term in brackets as the first term in a cumu- 
lant expansion. Re-exponentiate to get 



Z21 



K" 



(12) 



where 



Kf =K S - 6g( 



K s 



dr /r\3-2(K s +-if«) 



dr /r\ 
a \a/ 



n2(K s 



rfr 



3-2(A- s +A- a + 5 i-) 



■ (13) 



To obtain the RG equations we split the integrals into 
two parts 



/r\P f ba dr_ /ry f°° dr /r 

Ja a Ja a ^ Jba a ^ a 



(14) 



where In b <C 1. Perform the integral in the 1st term 
on the right hand side; rescale variables in the 2nd term 
so that the range of integration runs from a to oo. This 
gives 



-off _ r> c~2 / K s 



Kf =K S - 6g( 



-3g 4 



K. 



dr / T\3-2(K S +K a ) 

AirUsJ J a \a/ 

2 - dr /r\3-2(K s +K a +^-) 



V 47TU< 



dr fry 
a \a/ 



(15) 



where we have introduced 



K s = K s - 6gj 



4iru 

~ gi = gi 6 MM) 

.94 = .94 tf-^M^h) 



(b ~ 1) - 3<tf 



47TM< 



(6-1) 
(16) 



Since In b <C 1, we can write 1 — 6 = In 6 and b x = 
1 + x In b. Then write K s = K s + d K s , g\ = g\ + d gi , 
gi = 54 + d gi, and define dl = In b. This gives us the 
RG equations 



dK s 
dl 

d gi 

rf_94 
d2 



-6tf 



if. 



3^4 



- [2-{K s + K a )] 9l , 



34- 



(17) 
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